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ABSTRACT 

We present results from an analysis of stellar population parameters for 7132 galaxies in 
the 6dFGS Fundamental Plane (FP) sample. We bin the galaxies along the axes, v\, v%, and 
^3, of the tri-variate Gaussian to which we have fit the galaxy distribution in effective radius, 
surface brightness, and central velocity dispersion (FP space), and compute median values of 
stellar age, [Fe/H], [Z/H], and [a/Fe]. We determine the directions of the vectors in FP space 
along which each of the binned stellar population parameters vary most strongly. In contrast 
to previous work, we find stellar population trends not just with velocity dispersion and FP 
residual, but with radius and surface brightness as well. The most remarkable finding is that 
the stellar population parameters vary through the plane (v\ direction) and across the plane 
(^3 direction), but show no variation at all along the plane (i>2 direction). The t>2 direction 
in FP space roughly corresponds to 'luminosity density'. We interpret a galaxy's position 
along this vector as being closely tied to its merger history, such that early-type galaxies with 
lower luminosity density are more likely to have undergone major mergers. This conclusion 
is reinforced by an examination of the simulations of Kobayashi (2005), which show clear 
trends of merger history with i>2. 

Key words: galaxies:fundamental parameters, galaxies: elliptical and lenticular, galaxies: 
evolution, galaxies: structure 



1 INTRODUCTION 

Early-type galaxies are known to lie along a plane in the 3- 
dimensional (3D) parameter space whose axes are r=log(7? e ), 
s=log(<r), and i=log(7 e ), where R e , a, and I e represent effective 
radius, central velocity dispersion, and effective surface bright- 
ness respectively. This i s commonly referred to as the Fundamen- 
tal Pl ane (hereafter FP; [Dressier et all 1 19871 : iDiorgovski & Davis] 
[l987h . The plane can be expressed in the form 



r — as + bi + c (1) 
where a, b, and c are observationally derived constants. In the case 



where all early-type galaxies follow the virial theorem and have a 
constant mass-to-light ratio, a — 2 and b = — 1. In contrast to this 
theoretical expectation, the observed values are found to be in the 
range 1.2 < a < 1.6 and —0.90 < b < —0. 75 across a wide range 



of op ti cal and near infrared wavebands (e . g.jLucev . Bower & Ellis 
199ll IPahre. de Carvalho & Diorgovskl 1 19981. ' 



iHvde & Bernardi 



20091 ITa Barbera et alj2010af) . Thi s contrast may be due a t least in 



part to stellar population variations JCappellari et al ]2006h . though 
it has bee n argued that nonhomology mu st also contribute to the tilt 
(see, e.g.. lTruiillo. Burkert & Bell2004) . 

Stellar population variations lead to changes in the mass-to- 
light ratio, and may do so in ways that are correlated with FP pa- 
rameters, leading to tilts of the FP. Such variations may also in- 
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troduce additional scatter in the relation. Understanding these cor- 
relations and scatter may illuminate the origins of the FP and the 
formation of early- type galaxies. It may also lead to a means of 
improving the accuracy and precision of the FP distance estimator. 

Several authors have investigated correlations between 

FP param eters and st ellar po pulation (here after SP) pa- 

and 



rameters. 



Nelan et alj 



ISmith. Lucev & Hudsonl 
lations between a and 



iThomas et all J20051) . 
all found strong corre- 
several different SP parameters. 
iTerlevich & Forbesl 1 20021) found that [Mg/ F e] increases with lu- 
minos ity. iForbes. Ponman & Browr] J 19981) . iReda. Forbes & Hail 
d2005h . and iGargiulo et al .1 ( 120091) . among others, found cor- 
relations between SP p arameters and residuals from the FP. 
lLa Barbera et al. U2010bh found correlations between SP gradients 
and a, stellar mass, and dynamical mass. SP trends with radius 
and surface brightness individually, however, remained largely 
unexamined until recently. 

In their series of four papers l Graves. Faber & Schiavonl 

l2009al lbl: iGraves & Fabe3l20ld : IGraves. Faber & Schiavonll2010l) . 
Graves, Faber, & Schiavon took the analysis of SP trends in FP 
space a step further than the earlier studies, by considering SP 
trends along all three dimensions of FP space. They investigated 
SP parameter trends with radius, velocity dispersion, and surface 
brightness residuals from th e FP for galaxies in the Sloan D igital 
Sky Survey Data Release 6 dAdelman-McCarfhv et alj|2008l ; here- 
after SDSS), and found clear trends of several SP parameters with 
velocity dispersion, but much weaker trends with radius and surface 
brightness. 

The authors found that while each of the SP parameters they 
studied increases with increasing velocity dispersion, there is lit- 
tle COTjelatiOT_bj3twj^en^ and radius. In their Paper 
II ( IGraves. Faber & Sch iavon 2009bj; hereafter GFS), they hypoth- 
esized that galaxies with similar physical propert ies but different 
merg er histories may vary widely in radius (e.g.. lRobertson et al .1 
120061) . The fact that the SP parameters are insensitive to radius 
can thus be understood as an indication that SP parameters are 
determined independently of merger history. One drawback of the 
Graves, Faber, & Schiavon analysis, however, is that the authors 
bin galaxies along directions that are not orthogonal in FP space, 
which may potentially cause or hide correlations between SP and 
FP parameters. 

In this paper, we investigate SP trends in FP space using data 
from the Six-degree Field Galaxy Survey (6dFGS). 6dFGS is a 
near infrared and optically selected dual redshift/peculiar veloc- 
ity survey, with redshifts for over 125,0 00 galaxies in the south- 
ern hemisphere djones etal.ll2004l2009l) . 9572 of those galaxies 
are included in our 'velocity sample' (hereafter, 6dFGSv), which is 
described in Campbell et al. (in prep.). For each of the galaxies in 
6dFGSv, we have redshifts and velocity dispersions derived from 
6dFGS, along with J-, H-, and K-band radii and surface bright- 
nesses derived from 2MASS photometry. We have fit the FP to this 
sample as descri bed by Magoulas, Springob, & Colless et al. (sub- 
mitted, hereafter lMSQ) . 

Our ultimate aim is to derive distances and peculiar veloci- 
ties for these galaxies, which will be used to characterize the local 
galaxy peculiar velocity field and provide constraints on cosmo- 
logical models. However, we would first like to explore whether 
the FP relation's utility as a redshift independent distance indica- 
tor can be improved by accounting for SP variations in the plane. 
We would also like to gain a better understanding of both what SP 
trends mean for the star formation history in early-type galaxies, 



and what the distribution of galaxies in FP space means for both 
stellar and dynamical evolution. 

To this end, we have derived values of four SP param- 
eters (age, [Fe/H], [a/Fe], and [Z/H]) for 7132 of the galax- 
ies in 6dFGSv. As described in Proctor et al. (in prep.), these 
were derived using Lic k indices, following the procedure of 
iProctor & Sansoml ( l2002T) . In Section 2, we discuss the dataset, and 
briefly describe the fitting of the Fundamental Plane, as well as the 
derivation of the SP parameters. In Section 3, we present an anal- 
ysis of the SP trends in FP space in which we bin the galaxies by 
position in FP space, calculate the median value of each SP pa- 
rameter within each bin, and then calculate the gradient of the SP 
parameter variation, which gives us the vector in FP space along 
which the parameter varies. In Section 4, we discuss the physical 
interpretation of the SP trends in the context of the 3D Gaussian 
distribution of galaxies in FP space. 



2 DATA 

Complete details of the sample selection and da ta reduction are pre- 
sented in Campbell et al. (in prep.) and iMsd but we summarize 
the relevant points here. The 6dFGSv sample consists of all 6dFGS 
early-type galaxies that meet the following criteria: spectral signal- 
to-noise ratio greater than 5, heliocentric redshift ZhMo < 0.055, 
log of velocity dispersion s > 2.05 (in units of log[km/s]), and 
near infrared magnitude brighter th an m, = 13.65. As explained 
in both Campbell et al. (in prep.) and lMsd "early-type galaxies" in 
this context inclu des spi ral bulges in cases for which the bulge fills 
the 6dF fibre. As lMSd shows, spiral bulges follow essentially the 
same FP as early-type galaxies. While there is a 0.04 dex offset in 
the FP zeropoint between the sample of ellipticals and the sample 
of spiral bulges, there is actually no difference in the total thickness 
of the plane between the total sample and the sample of ellipticals. 
Including spiral bulges has no impact on the thickness of the plane. 

The apparent magnitudes used in this selection are taken from 
the Two Micron All-Sk y Survey (2MASS) Extended Source Cata- 
log djarrett etaL IHH). We have derived radii and surface bright- 
nesses for three different overlapping samples of 6dFGSv, corre- 
sponding to J-, H-, and K-band, with slightly different limiting 
magnitudes. Because the J-band sample offers photometric param- 
eters with the smallest errors, that is our preferred passband, and we 
have fit the FP in J-band. All photometric parameters used through- 
out this paper are in J-band. 

We have derived velocity dispersions for each of these galax- 
ies from their 6dFGS spectra. We have also derived half-light radii 
and surface brightnesses from their 2MASS J-band photometric 
images. Surface brightness is defined here as the average surface 
brightness interior to the half light radius. The angular radii have 
been converted to phy sical ra dii using the redshift distance to the 
galaxy, as explained in lMSd As mentioned in Section 1, we con- 
vert the physical radius, velocity dispersion, and surface brightness 
into logarithmic form, and write them as r, s, and i respectively. 

2872 of the galaxies are in groups or clusters, and we use the 
redshift distance of the group or cluster in such cases, where the 
group redshift is defined as the median redshift among all galaxies 
in the group. This is done because the galaxies within a group will 
tend to be at roughly the same distance, and the systemic redshift 
distance of the group offers a better estimate of the distance to the 
galaxy than the individual galaxy reds hift dis tance. Further details 
on the grouping algorithm are found in lMSd 

The initial morphological selection, described by Campbell et 
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al. (in prep.), involves matching the galaxy spectrum to a sample of 
spectral templates, and retaining only galaxies whose spectra match 
those of early-type galaxies. This selection allows for the inclusion 
of spirals, but is lik ely to do so only in cases for which the bulge 
fills the fibre. |MSC] describes how we subsequently inspected im- 
ages of each of the galaxies by eye, and eliminated cases with ob- 
vious problems, such as irregular morphologies. Spirals were only 
eliminated in cases for which either the fibre included some con- 
tribution from the disk or the galaxy was edge on and included a 
visible dust lane. 



2.1 Fitting the Fundamental Plane 

iMSd describes how we fit the FP using a maxim um likelihood 
method that closely follo ws the procedure of EFAR dColless et al .1 
200ll: Saglia et al. 200 1|). The procedure is explained in detail in 
Colless et alj feOOll) Section 3. It involves fitting the observed 
structural parameters to a 3D Gaussian in FP space. That is, we 
assume that, when plotted in r — s — i space, the galaxies follow 
a 3D Gaussian distribution, for which the two longest axes of the 
3D Gaussian define the Fundamental Plane, and the shortest axis of 
the 3D Gaussian is orthogonal to the plane. By construction, one of 
these axes has both an r and an i component, but no s component. 
The other two axes have components in all three dimensions. 

Given this 3D Gaussian functional form, the probability den- 
sity for the ith gala xy, P(xi), can then be computed according to 
IColless et al] d200lh Equation 1, where is the ith galaxy's po- 
sition in FP space relative to the axes of the 3D Gaussian. Given 
each galaxy's position in FP space, we then fit the orientation of 
the 3D Gaussian's axes (and thus the zeropoint and slopes a and 
b of the FP), by finding the orientation of the Gaussian that max- 
imizes the product of P(xi)'s for every g alaxy in the sample (the 
'likelihood'). (See Equations 2 and 6 in IColless et alj|200lL and 
the corresponding explanation in Section 3 of the paper.) The ac- 
tual orientation of the FP that gives the maximum likelihood is 
found by searching the multi-dimensional parameter space with 
a non-derivative multi-dimensional optimization algorithm called 
BOB YQA ( Bound Optimization BY Quadratic Approximation; 
|Powelll2006l) . 

The assumption of a 3D Gaussian distribution is motivated on 
purely empirical grounds. There is no obvious theoretical reason 
for one to expect that galaxies would follow such a distribution 
in FP space. However, when we fit to this model, the total likeli- 
hood of the fit (that is, the product of probability densities for every 
galaxy in the sample) is indistinguishable from that of mock cata- 
logs that were generated under the ass umptio n of 3D Gaussianity, 
suggesting a good fit. Additionally, as IMSd shows, the distribu- 
tions of individual parameters closely matches those generated by 
the mock catalogs as well. One might infer that the 3D Gaussian 
model implies that the luminosity function peaks at some interme- 
diate luminosity, and symmetrically falls off at fainter luminosities. 
This is not the case, however, as our selection cuts slice through FP 
space close to the center of the 3D Gaussian. Thus, dwarf galaxies 
are largely excluded, and the shape of the distribution of galaxies 
at the fainter end of the FP is unobserved. 

The best fit coefficients to the J-band FP r — as + bi + c 
are a = 1.524 ± 0.026, b = -0.885 ± 0.008, c = -0.329 ± 
0.054, where r, s, and i are in units of log[kpc/h], log[km/s], and 
log[Ls„ n /pc 2 ] respectively. (Note: The 'h' in kpc/h refers to the 
Hubble constant, in units of 100 km s _1 Mpc -1 . For the purpose 
of angular unit conversion, a flat cosmology of fi m = 0.3 and 



Axis of 3D Gaussian 
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i 




0.494 


-0.752 


0.437 


"2 


0.663 


0.000 


-0.749 


"3 


0.563 


0.659 


0.498 



Table 1. Transformation matrix between V\ — v% — V3 and r — s — i. 



Oa = 0.7 is assumed, though the specifics of the assumed cosmol- 
ogy affect the FP fit very weakly.) 

Following IColless et alj fcOOlh . we refer to the three axes of 
the 3D Gaussian as Vi, V2, and V3. The unit vectors along these axes 
are related to r, s, and i by the FP slopes a and b (from Equation 1 
of this paper) as follows: 

vi = f — as — bi, 

v 2 = f + l/b, (2) 
v 3 = -r/b- (l + b 2 )s/(ab) + i 

This closely follows lColless et al] j200lh . though there is one 
small difference, as we define surface brightness in log[Ls„ n /pc 2 ] 
units rather than magnitude units. Given our measured values of 
a — 1.524 and 6 = —0.885, we then have the transformation 
matrix between r — s — i space and V\ — v-z — vz space, which we 
provide in Table 1. That is, we provide the r-, s-, and i-components 
of unit vectors along vi, 1)2, and 773 and vice versa, (e.g., the unit 
vector along the Vi direction has length 0.494 in the r-directions, 
while the r unit vector correspondingly has length 0.494 in the in- 
direction.) vi is the shortest axis of the 3D Gaussian, orthogonal 
to the plane. It increases with increasing r and i, but decreasing s. 
V2, the longest axis of the 3D Gaussian, increases with increasing 
r and decreasing i, but has no s-component. W3 is the shorter of the 
two axes within the plane, and it increases with increasing r, s, and 
i. 

As explained by IMSd the longest axis of the 3D Gaussian 
(172) has no s-component by constru ction. We did perform one nine 
parameter fit (described in detail bv lMSd) in which we allowed V2 
to include an s-component. The best fit value of the s-component 
of the V2 unit vector is then -0.080. And so, in the nine-parameter 
fit scenario, the central value of the transformation matrix in Table 
1 becomes -0.080, and the other values in the matrix shift slightly 
so that every row and column is normalized to unity. Since this 
represents a very small shift, we exclude the s-component of the V2 
vector for all other fits. 

We also illustrate the directions of the vi, V2, and 173 axes rela- 
tive to the axes of physical parameters r, s, and i in Figure lQ. Also 
shown are the mass (M), luminosity (L), and mass-to-light ratio 
(M/L) directions, which we expl ain in more detail in Section 3.3. 

Following the convention of lMSd we refer to vi, 1)2, and V3 
as "through the plane", "along the plane", and "across the plane" 
respectively. 



1 This plot is an interactive 3D visuali zation, undertaken with custom C- 
code and the S2PLOT graphics library (games et al. 2006|). Interactive 3D 
figures, which can be accessed by viewing the version of this paper found 
in the ancillary files on astro-ph with Adobe Reader Version 8.0 or higher, 
were created using the approach described in iBarnes & Fluke] 120081) . 
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Figure 1. A 3D representation of the directions of the axes of the 3D Gaus- 
sian to which the FP was fit, v%, V2, and 113. The v%, V2, and ^3 axes are 
given in blue, with the axis vectors drawn proportionally to the length of the 
Gaussian (<ri, 0-2, and (73) along these three axes. The wireframe ellipses 
also illustrate the length of the Gaussian along the three axes. We note that 
vi is the shortest of the axes, and is thus orthogonal to the FP, while V2 
and V3 are directions within the FP. Also shown in red are the "mass", "lu- 
minosity", and "mass/luminosity" directions, as described in Section 3.3. 
Readers using version 8.0 or higher of Adobe Reader can enable interac- 
tive, 3D views of this plot by mouse clicking on the version of this figure 
found in the ancillary files. Once enabled, 3D mode allows the reader to 
rotate and zoom the view using the mouse. 



2.2 Derivation of stellar population parameters 

The x 2 -fitting procedure of IProctor & Sanso ml j2002l) was used to 
measure the derived parameters: log(age), [Fe/H], [Z/H] and [E/Fe] 
(which we hereafter refer to as [a/Fe], or the 'a' abundance ra- 
tio; see lThomas. Maraston & Ben der 2003 for details). Briefly, the 
technique for deriving these parameters involves the simultaneous 
comparison of as many observed indic es as possible to the model 
single stellar populations (SSPs) of iKorn. Maraston & Thomas! 
(2005). The best fit is found by minimising the square of the de- 
viations between observations and models in terms of the observa- 
tional errors (i.e. \ 2 )- 

The rationale behind this approach is that, while all indices 
show some degeneracy with respect to each of the derived parame- 
ters, each index does contain some information regarding each pa- 
rameter. In addition, such an approach should be relatively robust 
with respect to many problems that are commonly experienced in 
the measurement of spectral indices and their errors. These include 
poor flux calibration, poor sky subtraction, poorly constrained ve- 
locity dispersions, poor calibration to the Lick system and emission 
line contamination. This robustness is of particular importance in 
the analysis of large numbers of pipeline-reduced spectra such as 
those of the 6dFGS which cannot be accurately flux calibrated and 
so are not fully calibrated to the Lick system. The method is also 
relatively robust with respect to the uncertainties in the SSP models 
used in the interpretation of the measured indices; e.g. the second 
parameter effect in horizontal branch morphologies and the uncer- 
tainties associated with the Asymptotic Giant Branch. It was shown 



Figure 2. Histograms of all four SP parameters. 



in lProctor. Forbes & Beaslevl J2004I) and lProctor et alj | |2005|) that 
the results derived using the \ 2 technique are, indeed, significantly 
more reliable than those based on only a few indices. 

Fitting was carried out using an iterative clipping procedure. 
Initially the data were fit and a 5a clip was applied. The data were 
then re-fit and a 3a clip was applied. The fitting and 3a clipping 
were then iterated until no 3a outliers were found. Errors in the de- 
rived parame ters were estimated using the Monte Carlo technique 
described in IProctor et al' I d2008h - We note that, as described in 
IProctor & Sansoml \20o5) . the relationship between [Fe/H], [Z/H] 
and [a/Fe] is fixed such that [Z/H]=[Fe/H]+0.942[a/Fe]. 

A quality parameter was also defined as the sum of the inte- 
gerised reduced-^ 2 and the number of clipped indices. Only data 
from galaxies with S/N per angstrom greater than 9 and quality pa- 
rameter of 10 or lower are used in our analysis. As a result of the 
above procedures we measured the stellar population parameters in 
7132 galaxies, each utilising 10 or more indices. 

The distribution of each of the SP parameters is given in Fig- 
ure 2. The distribution of measurement errors on each of the SP 
parameters is given in Figure 3. The values for each parameter for 
each individual galaxy will be presented in Proctor et al. (in prep.) 



3 VARIATION OF STELLAR POPULATION 
PARAMETERS ACROSS AND THROUGH THE 
FUNDAMENTAL PLANE 

3.1 Global trends with physical parameters 

To examine the trends of SP parameters in FP space, we first bin 
the galaxies along the vi — — vs coordinate axes. That is, we 
set up bins in FP space with boundaries along those axes, of width 
0.1 in the vi direction, 0.2 in the V2 direction, and 0.2 in the W3 
direction. The bins are narrower along v± because we wish to have 
a comparable number of bins along each dimension. In each bin, 
we calculate the median value of each of the SP parameters. After 
removing all bins with fewer than 5 galaxies, we are left with 92 
bins. 

We now consider two approaches to assessing trends in FP 
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Figure 3. Histograms of the statistical uncertainties on all four SP parame- 
ters. 



space, using the median values in each bin. The first approach in- 
volves plotting each individual SP parameter against each FP pa- 
rameter. The problem with this type of analysis is that r, s, and i 
are each correlated with one another, and it is difficult to determine 
to what extent a trend with a particular FP parameter is merely an 
artifact of that parameter's trend with another FP parameter. 

Nevertheless, we have plotted the global variation of each SP 
parameter with each individual FP parameter in Figure 4. In this fig- 
ure, we show the median values of each of the four SP parameters 
vs. the corresponding r, s, and i values of each bin for our dataset. 
The figure includes a least squares fit to a linear trend for each of 
the individual SP-FP trends. We also plot a dashed line represent- 
ing a fit to a combination of "directional derivatives", that will be 
explained in Section 3.3. The figure also shows the R 2 correlation 
coefficient for each plot. All four SP parameters are seen to have 
a stronger correlation with velocity dispersion than with radius or 
surface brightness. 

This figure can be compared with Figures 4-6 of lGFS[ (N ote: 
We refer to [Mg/Fe] as [o/Fe]. We also use [Z/H], while iGFSl use 
[Mg/H], but the two quan tities are nearly identical.) Our results 
show agreement with lGFSl and other authors, in finding a positive 
correl ation between velocity dispersion and each of the SP param- 
eters. IGFSl find a weak positive correlation between r and both 
[Fe/H] and [Mg/H]. We find a similar positive correlation among 
these parameters (with [Z/H] in place of [Mg/H]), which may be 
slightly more pronounced in our data than in the GFS data. Addi- 
tionally, IgfsI claim a mild correlation between surface brightness 
and some of the SP parameters, but there are no such trends appar- 
ent in our Fig. 4. 



3.2 Variations along the axes of the 3D Gaussian 

As mentioned in the previous section, the correlations between r, 
s, and i complicate the interpretation of Figure 4. A more useful 
approach for displaying the data is to plot the full 3D distribution. 



Figure 4. Results of the stellar population modeling by bin in FP space. 
Each point is the median value of one of the SP parameters in a bin in FP 
space, plotted against the corresponding r, s, or i value at the center of that 
bin. We also superimpose the best lit regression line to the plotted points 
(solid line) and a best fit line to a set of directional derivatives for r, s, and i 
(dashed line, see Section 3.3). The R 2 correlation coeff icient is given in the 
upper left corner of each plot. As (Proctor et al. 2008) shows, each of the 
possible SP values is quantized, though the quantization is most extreme for 
[o/Fe], 



In Figures 5, 6, 7, and fff], we show 3-dimensional FP space varia- 
tion of A —log(age), [Fe/H], [«/Fe], and [Z/H] respectively. When 
plotted this way, one can see more complex trends than those ob- 
served in Fig. 4. 

The second approach we take to assess the SP trends allows us 
to examine this complexity: We examine the trend of each SP pa- 
rameter with each FP parameter when the other two FP parameters 
are held constant. This approach makes use of the partial deriva- 
tives dS/dJ-, where S is used as a shorthand representation for 
any of the SP parameters and T is used as a shorthand representa- 
tion for any of the FP parameters. 

We calculate these partial derivatives in the following way. 
First, we assume that the variation of each of the SP parameters 
in FP space can be fit by a least squares fit to a straight line. We 
perform such a linear regression along the Vi, V2, and v$ direc- 
tions individually, using the median values of each SP parameter 
in every bin for each fit. For example, in fitting the variation of 
A —log(age) along the Wi direction, we fit the partial derivative of 
log(age) with Vi according to: 



dA 



NHviiAi — EwiiEAi 



(3) 



- (£vi«) 2 

where Vu and Ai are, respectively, the v\ position of the ith bin 
and the median value of A of the galaxies in the ith bin. The 
summation is performed over all N — 92 bins containing 5 or 
more galaxies. By then making the corresponding calculations for 



These are interactive 3D figures, generated in the same manner as Figure 



1. 
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Figure 5. Variation of log(age) across the Fundamental Plane, in 3D. Each 
sphere represents a bin in FP space, including 5 or more galaxies. The 
sphere is placed at the midpoint of the bin's r — s — i coordinates, color- 
coded so that redder colors represent older ages, and bluer colors represent 
younger ages, as given by the color scale on the right of the plot. The size 
of the sphere scales with the logarithm of the number of galaxies in the bin, 
as given by the scale established by the black spheres on the side of the 
plot. The number labeling each of the black spheres is the logarithm of the 
number of galaxies in a bin represented by a sphere of that size. As with 
Figure 1, readers using version 8.0 or higher of Adobe Reader can enable 
3D interactive views of this plot by mouse clicking on the figure. 



[Fe/H] 




Figure 6. Same as Figure 5, but with [Fe/H] rather than log(age). Redder 
colors indicate higher values of [Fe/H]. 

the V2 and v$ directions, we produce a vector of age variation in 
Vi — V2 — V3 space. 

While the components of this vector are partial derivatives, the 
vector itself can be thought of as the gradient of age in vi — V2 — «3 
space, \jA. In the next subsection, we generalize this approach 
to include the components of the gradient along directions other 
than Vi,V2, and v$. In such cases, it no longer makes sense to de- 
scribe the components of the gradient as merely 'partial deriva- 
tives' . Rather, we will refer to them as 'directional derivatives' . The 
directional derivative of A with respect to vi, which we write as 




Figure 7. Same as Figure 5, but with [ct/Fe] rather than log(age). Redder 
colors indicate higher values of [ct/Fe]. 



[Z/H] 




Figure 8. Same as Figure 5, but with [Z/H] rather than log(age). Redder 
colors indicate higher values of [Z/H]. 



SJviA, is the change in A along the vi direction, per unit change 
in Likewise, \JfA is the change in A per unit r, and X/rhA the 
change in A per unit m — \og(mass). 

We have computed directional derivatives for each of the SP 
parameters: log(age), [Fe/H], [a/Fe], and [Z/H], which can be 
found in Table 2. We also provide the statistical errors on each of 
the directional derivatives, e, as well as the absolute value of the ra- 
tio of each directional derivative to its own error, \ (e.g., for v -pA, 
X = | V;r Trends with significance \ > 3 are bolded. 

In using this method, we have implicitly assumed that, for any 
given SP parameter, there is a direction in FP space along which 
that parameter increases linearly. To test this hypothesis, we have 
also fit a quadratic curve to each of the SP trends along vi, i>2, and 
V3 . In every case, we find that the quadratic coefficient is consistent 
with zero, to within the statistical errors. The assumption of linear 
variation in FP space thus seems justified. In fact, even if there were 
minor deviations from linearity, our method would still be sufficient 
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-1.47 


0.12 


12.25 


0.37 
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0.05 


4.80 


0.07 


0.13 


0.54 


''2 


-0.04 


0.04 


1.00 


0.05 


0.02 


2.50 


-0.01 


0.01 


1.00 


0.05 


0.03 


1.67 


V3 


0.08 


0.09 


0.89 


0.26 


0.04 


6.50 


0.16 


0.02 


8.00 


0.46 


0.04 


11.50 


r 


-0.70 


0.08 


8.75 


0.37 


0.06 


6.17 


-0.03 


0.03 


1.00 


0.32 


0.07 


4.57 


s 


1.16 


0.11 


10.55 


-0.11 


0.08 


1.38 


0.29 


0.04 


7.25 


0.25 


0.10 


2.50 


i 


-0.57 


0.08 


7.13 


0.25 


0.05 


5.00 


-0.02 


0.03 


0.67 


0.22 


0.06 


3.67 


in 


0.32 


0.05 


6.92 


0.03 


0.03 


0.88 


0.11 


0.02 


6.44 


0.16 


0.04 


3.87 


I 


-0.39 


0.04 


11.01 


0.20 


0.03 


7.62 


-0.02 


0.01 


1.19 


0.17 


0.03 


5.65 


m. — I 


0.60 


0.04 


14.51 


-0.14 


0.03 


4.72 


0.11 


0.02 


6.96 


-0.01 


0.04 


0.18 



Table 2. Stellar population trends in FP space. 



to illuminate these qualitative relationships between the SP and FP 
parameters. 

One remarkable feature of the results in Tables 2 is the lack 
of variation of any of the SP parameters along the V2 direction, 
the long dimension of the FP. All of the SP parameters vary along 
a direction that is a superposition of the vi ('through the plane') 
and V3 ('across the plane') axes. Age variation is almost entirely 
through the plane, while [Z/H] variation is almost entirely across 
the plane, and [a/Fe] and [Fe/H] are superpositions of the two. The 
lack of variation of any of the SP parameters along i>2 is a major 
result of this paper, and is examined in more detail in Section 4. 

We should also note here that v± is a quantity that has been 
studied by other authors. It is simply the residual from the plane, 
measured along the direction orthogonal to the plane. Several au- 
thors have examined correlations between various parameters and 
FP residual. In some cases, this is measured as the residual along 
a different dimension, such as radius. However, this scales with vi 
to within a constant scale factor, so increasing v± is proportional 
to increasing residual in r. lGargiulo et al.l j2009h . for example, find 
an anticorrelation between age and residual in r, as well as an an- 
ticorrelation between [a/Fe] and residual in r. This is consistent 
with our finding of an anticorrelatio n between both of thes e quan- 
tities and v±. However, in contrast to lGargiulo et alj we do 
not find that the anticorrelation with [a/Fe] is stronger than the one 
with age. 



3.3 Variations with physical parameters 

As mentioned in Section 3.2, we wish to generalize the derivation 
of the directional derivative to directions other than Vi,V2, and vs. 
We would like to examine how the SP parameter variations in vi — 
V2 — V'a space translate to variations in r — s — i space. In the case 
of the directional derivative with v\, S/eiA is exactly the same as 
dA/dvi, as given by Equation 3. While in the case of mass, the 
expression for the directional derivative of age with mass, SJmA, is 
a linear combination of partial derivatives with vi, V2, and 113. This 
is also the case for the directional derivative with luminosity, mass- 
to-light ratio, radius, velocity dispersion, and surface brightness. 

Table 1 gives the Vi, V2, and V3 components of r, s, and i. 
This provides the coordinate transformation from \/$iS, V«2<5, 
and V*3<S to S7fS, \JgS, and S/iS. We are also interested in calcu- 
lating how each of the SP parameters varies with dynamical mass, 
luminosity, and mass-to-light ratio. If we assume homology, then: 



I = 2r + i + c 2 



(5) 



where m = log(mass), I = \og(luminosity) , and ci and C2 are 
normalization constants. Subtracting these equations, we express 
the logarithm of the mass-to-light ratio as 



m — I = —r + 2s — i + ci — C2 



(6) 



We then wish to derive \JmS, Vf^>> and \/ m ^ t S. As previously 
mentioned, we must be careful about how we define the directional 
derivative for a quantity such as mass, which does not represent any 
of the basis vectors in FP space. We are defining \/fhS to mean 
the change in an SP parameter <S per unit m along the gradient 
of m in FP space, \jm. This direction is xjm = r + 2s, or the 
direction along which, for every increase (Sr) in r of one unit, there 
is a corresponding increase (Ss) in s of two units. To normalize 
this vector to a change (Sm) in m of 1 , we should actually divide 
by 5, because we require Sm — Sr + 28s, and 8s — 28r. Thus, 
Sm — Sr + 454r — 5Sr, and Sr — 1/5. 

Thus, the directional derivative of the SP parameter 5, along 
the normalized gradient of m is: 



_ 1 dS 2dS 
5 or 5 as 



(7) 



This is the change in 5 per unit change in m, provided that m 
is changing along its gradient in r — s — i space (along the direction 
(5r,5s,Si) = (+l,+2,0)). 

We similarly derive: 



5 dr 5 di 



Vm - 1 ~ 6 dr + 3 ds ~ 6 di 



(8) 



(9) 



r + 2s + ci 



(4) 



The resulting relationships between each SP parameter and each 
structural parameter are given in Table 2. We also include statis- 
tical uncertainties e v g and the ratio between and e v s. We 
note that the m, I, and m — I directions in FP space that we have de- 
rived here correspond to the directions shown in Figure 1 (written 
as log(A/), log(L), and log(M/L) respectively). 

We have also taken the individual directional derivatives 
Vs>5, and Vi<S> an d computed the inferred variation of the SP pa- 
rameters for a set of bins matching our bins' positions in FP space, 
then fit a regression line to those points. The best fit lines are shown 
as the dashed lines in Figure 4. The extremely close match between 
these best fit lines and the solid lines (from fits that do not assume 
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there is a single direction in FP space along which the SP parameter 
varies linearly) offers further evidence that our linear fits are a good 
match to the real SP parameter variation. 

As Table 2 shows, there are substantial differences between 
the different SP parameters in terms of how they vary with the 
FP parameters. Three of the four SP parameters increase with in- 
creasing velocity dispersion. The outlying case is [Fe/H], which 
shows no dependence on s for fixed r and i, despite the fact that, 
as discussed earlier, there is a global correlation between [Fe/H] 
and s when r and i are allowed to vary. Dependence on r and i 
varies from parameter to parameter, with [Z/H] and [Fe/H] increas- 
ing with increasing r and i, age increasing with decreasing r and i, 
and [a/Fe] independent of r and i. Age varies most strongly with 
Wi, [Fe/H] most strongly with vi and r, [a/Fe] most strongly with 
s, and [Z/H] most strongly with W3. 

We note that some of these trends have been identified by other 
authors as well. For example, the relati onship between age and FP 



residu al, which we call vi, was noted bv lForbes. Ponman & Browr 



1998). And as we noted in the previous subsection. lGargiulo et al 



2009) found relationships between FP residual and both age 



and [a/Fe] that are consistent with ours, at least in terms of 
the sense of the correlation. We also note that our estimate of 
dA/ d(m — I) is 0.60, which contrasts with ~ 0.75 (as estimated 
from lBruzual & Charlotl2003l. Figure 3), which is predicted by stel- 
lar population models. The difference may well reflect differences 
between dynamical mass and stellar mass. 

In summary, our results for the global trends of the SP param- 
eters with respect to the FP parameters agree in broad terms with 
the trends observed by other authors, in that each SP parameter is 
seen to be positively correlated with velocity dispersion. However, 
we have also taken the additional step of deriving the directional 
derivatives of each of the SP parameters with each of the FP pa- 
rameters (Table 2). That is, we have calculated the dependence of 
each SP parameter on each FP parameter by computing its gradient 
in FP space. When analyzed in this way, it is no longer the case that 
all four SP parameters depend more strongly on velocity dispersion 
than any other FP parameter. Most interestingly, the vectors along 
which the SP parameters vary are closely aligned with the axes of 
the 3D Gaussian that defines the FP, with age varying almost en- 
tirely with ui, [Z/H] varying almost entirely with v$, and none of 
the SP parameters varying along V2 . 

3.4 Comparing 'slices' of 6dFGS FP with those of SDSS 

We have shown how examining the full 3-dimensional distribution 
of stellar population variations in FP space gives one a clearer pic- 
ture of the SP trends than one would get from collapsing the trends 
down to a 2-dimensional distribution. We now present anot her set 
of plo ts that we can compare to the SDSS results presented in lGFSl 
IGFSI Figures 7-10 show the variation of the SP parameters within a 
given slice of the FP: below the plane, within the plane, and above 
the plane. We produce similar figures for our data in our Figures 
9-12. 

Figure 9 shows our FP bins distributed in three slices. The 
left panel shows the distribution of bins directly "below" the plane. 
That is, it shows the bins with central v i valu es of -0.1, equivalent 
to negative Ai in the nomenclature of iGFSl The two axes are r' 
and s', where r' is the direction within the FP along which r in- 
creases but s remains constant and s' is the direction within the FP 
along which s increases but r remains constant. The middle panel 
is the equivalent plot for the 'midplane', which includes the bins 
with central wi = 0, and the right panel is the equivalent plot for 



bins with central vi = 0.1, equivalent to positive Ai in the IGFSl 
nomenclature. Our binning scheme also includes a few bins with 
central vi = —0.2 and central vi = 0.2, but they include very few 
galaxies, and are not plotted here. 

In each plot, the area of each circle is proportional to the loga- 
rithm of the number of galaxies in the bin, and the color of the cir- 
cle scales with log(a<;e), with redder points corresponding to older 
galaxies. The main trend seen here is the large variation occurring 
between different slices of vi, with the bins 'above' the plane in- 
cluding more young galaxies. 

In Figure 10, we present the equivalent plot for [Fe/H]. In 
this case, [Fe/H] increases with increasing Vi. We also observe that 
[Fe/H] increases with increasing r' and s'. In Figure 1 1, we present 
the equivalent plot for [a/Fe]. [a/Fe] is positively correlated with s' 
at low vi, whereas [a/Fe] is more consistently low and only weakly 
dependent on s' at high Vi. In Figure 12, we present the equivalent 
plot for [Z/H]. [Z/H] increases sharply with increasing s', though 
the trend is somewhat muted for the high vi slice relative to the 
other two slices. 

Each of these trends closely tracks the SP trends that lGFSl ob- 
serve in the SDSS Fundamental Plane, when the data are plotted in 
the sa me manner (their Fig. 7-10). (We do note just two differences: 
IGFSl find a much stronger trend between age and s' and a weaker 
trend between [Fe/H] and r' than we do.) This suggests that the 
underlying trends in the SDSS and 6dFGS datasets are very nearly 
the same. 



4 DISCUSSION 

iGraves, Faber & Schiavor] J2010h characterize the SP parameter 
variation in FP space as 'the 2D family of early-type galaxy stellar 
populations and their structural properties'. This is contrasted with 
the ID mass sequence of gal axies, which was t he f ocus of earlier 
studie s. Earlier work, such as lNelan et al.l d2005h and lThomas et al.l 
d2005l) , examined the variation of the SP parameters with a, which 
was used as a proxy for galaxy mass. To first order, one can imagine 
these ID relationships as variations along the a projection within 
theFP. 

Graves, Faber, & Schiavon improve on this analysis by also 
exploring SP variations through the plane. The second dimension in 
the 2D family of early-type SPs represents residuals from the FP. SP 
variations along this dimension can be interpreted as evolutionary 
differences, with aging stellar populations evolving in mass-to-light 
ratio. Our results show, for example, that age vari es more strongly 
than th e other SP parameters with FP residual. As lGraves & Faberi 
j2010l) point out, however, the thickness of the plane may be better 
explained by genuine structural differences than by the fading of 
stellar populations. 

Table 3 of IGraves & Faber! l l2010h summarizes the qualitative 
relationships between each SP parameter and each FP parameter 
for their dataset. However, as explained in the previous section, the 
three directions in FP space that the authors consider are not or- 
thogonal. They are what we refer to as r' , s', and Ai, where r' 
and s' are as defined in the previous section (the direction within 
the plane along which r increases but s remains constant and the 
direction within the FP along which s increases but r remains con- 
stant, respectively) and Ai is the residual from the plane along the 
i direction. 

This scheme is similar to, but not quite the same as, our or- 
thogonal vi — V2 — Vz basis set. r' is exactly identical to vi. Ai is 
defined somewhat differently from vi , but the two are equivalent to 
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Figure 9. Each panel shows the distribution of bins (described in Section 3.1) across r' and s' for a different slice of the FP. The left 
panel is the vi < slice ('below the plane'), the middle panel is vi ~ ('the midplane'), and the right panel is the vi > slice 
('above the plane'), r' and s' are as described in Section 3.4: r' is the value of r within the plane for fixed s, and s' is the value of s 
within the plane for fixed r. The area of each circle is proportional to the logarithm of the number of galaxies in the bin. The color of 
each circle represents the median value of log(age) in the bin, as given by the color scale shown on the right. Redder colors correspond 
to older ages, and bluer colors correspond to younger ages. 
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Figure 10. Same as Figure 9, but for [Fe/H] rather than age. Redder colors correspond to higher values of [Fe/H]. 



one another to within a constant scaling factor. The one dimension 
that Graves et al. use that is different from any of our basis vectors 
is s' . It includes a significant v$ component, but also has a «2 com- 
ponent. However, since none of the SP parameters varies with V2, 
we would expect that any SP trend with vs would also be apparent 
in s' . 

And indeed, if one compares iGraves & Faberl feOlOh Table 3 
with Table 2 from this paper, one gets remarkable agreement, from 
substituting Ai for vi, r' for V2, s' for «3, [a/Fe] for [Mg/Fe], and 
[Z/H] for [Mg/H]. If one classifies any trend detected in our data 
at a significance of less than 3<r as 'null' in the nomenclature of 
Graves & Faber, then the only differences between the two tables 
upon making such substitutions are that we find no statistically sig- 
nificant trends between age and ^3, nor between [Z/H] and vi. So 
despite the differences in the binning scheme and despite the fact 
that Graves et al. stack the spectra in each bin, whereas we de- 
rive SP parameters for individual galaxies, we find broadly similar 
trends for each of the SP parameters. 

While the analysis by [Graves & Faberl (f2O10h focuses on the 
'2D' variation of SP parameters, we would actually like to focus in 
particular on the third dimension in FP space, the dimension along 
which none of the SP parameters seems to vary. As we show in 
Table 2, the vectors along which the SP parameters vary seem to 
be more closely aligned with the vi and 113 vectors than they do 
with the vectors of any simple physical quantity. Age varies almost 
purely along Wi, [Z/H] varies almost purely along v$, and [Fe/H] 



and [a/Fe] vary along superpositions of v\ and 113, with no compo- 
nent in V2 . 

This is a remarkable result, as we had no physical motivation 
in choosing the directions of these vectors, allowing the data to 
determine the axes of the fitted 3D Gaussian. However, because 
the SP parameters are so closely aligned with the axes of the 3D 
Gaussian, it now seems likely that there is in fact some physically 
meaningful reason for the V2 and v$ axes to be oriented in these 
particular directions in FP space. What we require is a hypothesis 
that explains both the distribution of galaxies in FP space, and why 
the SP parameters vary along the 3D Gaussian. 



4.1 V2 and merging history 

In Section 5 of lGFSl the authors discuss the question of why the SP 
parameters apparently vary with velocity dispersio n, but not radius. 
They suggest that, as seen in Af-body simulat ions jRobertson et al .1 
120061 : iBovlan-Kolchin. Ma & Ouataertl [20051) . galaxies of compa- 
rable mass and luminosity with different merger histories can have 
wildly different radii and surface brightnesses. Velocity dispersion, 
on the other hand, is relatively independent of merger history. The 
argument is that if SP parameters are independent of merger history, 
then it would follow that SP parameters would vary with velocity 
dispersion, but not radius. 

In this work, however, we find that there is in fact variation 
of SP parameters with radius, even within the plane, as most of 
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Figure 11. Same as Figure 9, but for [a/Fe] rather than age. Redder colors correspond to higher values of [ce/Fe]. 
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Figure 12. Same as Figure 9, but for [Z/H] rather than age. Redder colors correspond to higher values of [Z/H]. 
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the SP parameters vary with V3, which includes an r component. 
We believe ou r resu lts can be explained by a somewhat modified 
version of the iGFSl hypothesis: The variations in merger history 
add significant scatter to the correlations between radius and the 
SP parameters, but they do not eliminate them completely. 

This idea is supported by the simulations of lKobavashil J2005h . 
Figure 5 of that paper shows how variable merger histories have 
little impact on the scaling relations of a with other physical pa- 
rameters of the galaxy, while adding significant scatter to scaling 
relations with r. However, despite the significantly added scatter, a 
residual correlation between r and the other galaxy parameters re- 
mains. There is no reason why this should not also be the case with 
respect to correlations between r and the SP parameters. 

We then offer the following scenario: As one moves along the 
113 axis towards increasing r, s, and i, one finds galaxies of in- 
creasing mass, luminosity, and metallicity. (This is because larger 
r, s, and i necessarily implies larger mass and luminosity, as 
log(mass) = r + 2s and log(luminosity) = i + 2r.) This could 
simply be driven by the increasing total mass of the system, or 
the increasing mass of the dark matter halo that seeded the galaxy. 
However, as one moves along the V2 axis, one finds galaxies of in- 
creasing radius and decreasing surface brightness associated with 
variations in merger history. Indeed, if the FP were precisely the 
virial plane, then V2 would be precisely luminosity /radius (lu- 
minosity density). We conclude then that luminosity density is de- 
termined by, or at least heavily influenced by, merger history, which 
in turn has little to no impact on the galaxy's stellar population. 

This hypothesis is also cons istent with the results of 
iTruiillo. Ferreras & de la Rosal d201lh . who examine samples of el- 
liptical galaxies at both z ~ and z ~ 1, and find that evolution 
in size is independent of stellar age. The authors suggest that this 



argues in favor of size evolution since z ~ 1 being driven more by 
dry mergers than by a "puffing up" scenario dFan et al.l2008l [20 1 Ol : 
iDamianov et al.| [2009i) in which growth in galaxy radius is due to 
the expulsion of gas by AGN or supernova-driven winds. 

If our hypothesis is correct, then galaxies are spread along the 
V2 and V3 directions by physically unrelated processes, related to 
secular migration and the galaxy's initial conditions respectively. 
It is then something of a 'coincidence' that these two processes 
distribute the galaxies along orthogonal directions in FP space. Of 
course, we cannot rule out small deviations from orthogonality that 
our fitting method is unable to detect. 



4.2 Simulations of merger history variations across the 
Fundamental Plane 

One approach we can take to investigate this issue further is to ex- 
amine the distribution of galaxies in FP space, as seen in simu- 
lations, for which the merger history of each galaxy is known. To 
that en d, we have examined the simulations presented in lKobavashil 
d2005h . In that paper, the author simulates the formation and evo- 
lution of 128 galaxies using a smoothed particle hydrodynamics 
method and a special purpose computer GRAPE (GRAvity PipE). 
Because both the FP parameters and formation mechanisms of each 
of these simulated galaxies is known, we can compare the merger 
histories to the galaxy distribution in FP space, and determine what 
impact merger history has on a galaxy's position in FP space. We 
describe this investigation below. 

The details of the GRAPE-SPH code are presented in 
iKobavashil l2004h . but we briefly recount the main features of the 
simulation here: The initial conditions are 74 spherical regions with 
Cold Dark Matter initial fluctuations, which produces 83 ellipti- 
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Figure 13. Fundamental plane relation for galaxies in the lKobavashlfeOOSl) 
simulation. The colors show the merging histories of the galaxies, following 
the classification scheme explained in Section 3.1 of that paper: [El] mono- 
lithic (red), [E2] assembly (magenta), [E3] minor merger (green), [E4] ma- 
jor merger (cyan), [E5] multiple major mergers (blue), and [Dl-5] dwarfs 
(black crosses). The plot shows that while there is variation in merger his- 
tory along the plane, there is no apparent trend between merger history and 
scatter off of the plane. 



cal galaxies and 45 dwarf galaxies. As well as the kinematics of 
the dark matter and gas particles, the relevant baryon physics (i.e., 
radiative cooling, star formation, chemical enrichment, and super- 
nova feedback) are included. Computing the particle distribution 
from z ~ 25 to z — 0, the time evolution of the internal structures 
of galaxies are predicted. 

Physical parameters are measured for each of t hese simu- 
lated g alaxies, and the scaling relations are examined in lKobavashil 
(2005). As Figure 8 of that paper shows, if one plots th e resulting 
parameters in K-space l lBender. Burstein &Faberlll992h . it appears 
that divergent merging histories increase the scatter of the Funda- 
mental Plane. 

We note, however, that the K-space representation does not 
account for the tilt of the FP. In Figure 13, we thus produce our own 
plot of the "edge-on" view of the FP for these simulated galaxies in 
FP space. We have fit the FP for the simulated galaxies in question, 
and find that they follow an FP relation of r = 1.30s — 0.54i — 1.03. 
Deviations from the line in Figure 13 thus represent deviations from 
the FP. Galaxies with all of the realized merging histories appear to 
follow the same FP relation. 

In Figure 14, we plot s vs V2, where Wa = 0.475r — 0.880i, 
as appropriate for the FP relation for the simulated dataset in ques- 
tion. Dwarf galaxies are excluded, as they nearly all have s < 2.0, 
and thus do not have counterparts in the 6dFGSv sample. This plot 
shows that there are clear correlations between merger history and 
position in FP space. The ellipticals that formed via monolithic col- 
lapse or the assembly of subgalaxies are preferentially found to 
have large velocity dispersions, and low values of V2, which corre- 
sponds to small radii and large surface brightnesses. For any fixed 
value of s, there appears to be a relationship between the value of 
V2 and the merger history of the galaxy, such that galaxies with 
large radii and small surface brightnesses are likely to have been 
formed by one or more major mergers. 

In Section 4.1, we hypothesized that variations in merger his- 
tory elongate the FP along the V2 direction. Figure 14 would seem 
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Figure 14. V2 vs. s for the elliptical galaxies in the Kobayashi 1 2005) sim- 
ulation. The colors are defined the same as in Figure 13. We omit dwarf 
galaxies from the plot, as they would not be included in 6dFGSv anyway. 
The plot shows clear trends between merger history and both 1)2 and cr. 



to confirm that this is the case for the simulations, as we do indeed 
see that for fixed s, there is a dependence of 1)2 on merger history. 

There is, however, one difference between these simulations 
and the real data of 6dFGSv that complicates our interpretation, 
and should be noted. The value of b in these simulations is -0.54, 
very close to the -0.5 value that corresponds to constant luminosity 
along the V2 direction, and quite different from the virial expecta- 
tion of -1. In the simulations, when one compares different galaxies 
at the same velocity dispersion, but different values of 1)2, the lu- 
minosities are nearly the same. Thus, the variations in merger his- 
tory do not lead to variations in luminosity, when one controls for 
s. This does not hold in the real data however. Since the 6dFGSv 
sample has b = —0.885, the galaxies with larger V2 have not just 
larger radii, but larger luminosi ties. 

Likewise, as discussed in ImscI essentially all authors have 
found a value of b in the range — 0.9 < b < —0.7 , regardless of 
waveband. The lower right panel o f lKobavashil l l2005l) Figure 5 con- 
firms the mismatch in the slopes of r vs. i between the simulations 
and real data. The reason for this discrepancy is unclear. It is possi- 
ble that this is because the secondary starbursts are rarely induced 
by mergers in the simulations. If secondary starbursts occur, the 
metallicity should be larger, which may break the mass-metallicity 
relation of galaxies. This should be addressed in future simulations 
with higher resolutions. 

Since we do not find b = —0.5 in our dataset, there remains 
something of a problem with the hypothesis that V2 represents vari- 
ations in luminosity density from divergent merger histories. What, 
precisely, does it mean for two galaxies to be "the same except 
for merger history"? In the simulations, luminosity is held nearly 
constant along V2. In the real data, however, what is held constant 
along «2? Neither mass nor luminosity is fixed along V2, though 
they both vary far more slowly along V2 than along W3. One possi- 
bility is that galaxies along V2 (with fixed Vi and 1)3) had the same 
initial mass at an earlier epoch, but have accreted different amounts 
of additional material. Since the slope of V2 also does not match the 
virial expectation of b = —1, galaxies with larger V2 not only have 
larger radii, but slightly larger mass-to-light ratios. We would thus 
conclude that galaxies which underwent major mergers would, on 
average, have somewhat greater mass-to-light ratios. 
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